Robustness of s-wave Pairing in Electron-Overdoped Ai_j,Fe2-xSe2 
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Using self consistent mean field and functional renormalization group approaches we show that s- 
wave pairing symmetry is robust in the heavily electron-doped iron chalcogenides (K, Cs)Fe2-2;Se2. 
This is because in these materials the leading antiferromagnetic (AFM) exchange coupling is between 
next-nearest-neighbor (NNN) sites while the nearest neighbor (NN) magnetic exchange coupling is 
ferromagnetic (FM). This is diff'erent from the iron pnictides, where the NN magnetic exchange 
coupling is AFM and leads to strong competition between s-wave and d-wave pairing in the elec- 
tron overdoped region. Our finding of a robust s-wave pairing in (K, Cs)Fe2-a;Se2 differs from 
the d-wave pairing result obtained by other theories where non-local bare interaction terms and 
the NNN J2 term are underestimated. Detecting the pairing symmetry in (K, Cs)Fe2-a;Se2 may 
hence provide important insights regarding the mechanism of superconducting pairing in iron based 
superconductors. 



I. INTRODUCTION 

The recent discovery of a new family of iron-based su- 
perconductors A(K,Cs,Rb)j^Fe2_a;Se?i"— has initiated a 
new round of research in this field. Remarkably, this 
new family shows distinctly different properties from 
other pnictide families: the compounds are heavily elec- 
tron doped, but their superconducting transition tem- 
peratures are high, at more than 40 K. For compar- 
ison, such large Tc's can only be reached in the op- 
timally doped 122 iron pnictides^. Importantly, both 
angle-resolved photoemission spectroscopy (ARPES)^"— 
and LDA calculations^"— show the presence of only elec- 
tron Fermi pockets located at the M point of the folded 
Brillouin zone (BZ). (Some signature of possible density 
of states at the F point is still under current debate; in 
any case this pocket, if present, is assumed to be very flat 
and shallow). ARPES experiments have also reported 
large isotropic superconducting gaps at these pockets^"—. 
The absence of hole pockets around the F point of the 
BZ provides a new arena of Fermi surface topology to in- 
vestigate the pairing symmetries and mechanisms of su- 
perconductivity proposed for iron-based superconductors 
from a variety of different approachesii"— . 

So far, the majority of theories for the pairing sym- 
metry of iron-based superconductors are based on weak 
coupling approachesi^"— i^"— i^i^i. Although there are 
discrepancies, the theories based on these approaches 
have reached a broad consensus regarding the pairing 
symmetries in iron-based superconductors: for optimally 
hole doped iron-pnictides, for example, Bao.6Ko.4Fe2As2, 
an extended s-wave pairing symmetry, called s^, is 
favoredi^ (the sign of the order parameter changes be- 
tween hole and electron pockets as potentially detectable 
through neutron scattering^^) , as a result of repulsive in- 
terband interactions and nesting between the hole and 
electron pockets. For extremely hole-doped materials, 
such as KFe2As2, the absence of electron pockets can 
lead to a d-wave pairing symmetryi^ with a low transi- 



tion temperature; for electron doped materials such as 
Ba2Fc2-2;Co3;As2, the anisotropy of the superconducting 
gap around the electron pockets in the state grows for 
larger electron doping and eventually the SC gap devel- 
ops nodes around the electron pockets due to the weak- 
ening of the nesting condition and the increase of dxy 
orbital weight at the electron pocket Fermi surfaces^^iiS. 
Finally, in the limit of the heavily electron doped case 
when the hole pockets vanish and only the electron pock- 
ets are left, the d-wave pairing symmetry may be favored 
g^gg^jjji8j4i|42_ rpj^g j^p^ chalcogenide AyFe2-a;Se2 belongs 
to the latter category and many theories based on weak 
coupling approaches have suggested that the pairing sym- 
metry should be d-wave as possibly detectable through 
characteristic impurity scattering^"—. 

A complementary approach based on strong coupling 
likewise predicts an s-wave pairing symmetry in the iron 
pnictides. Two of us showed that the pairing symme- 
try is determined mainly by the next-nearest-neighbor 
(NNN) AFM exchange coupling J2 together with a renor- 
malized narrow band widt h^^'^^ . The superconducting 
gap is close to a cos kx cos ky form in momentum space 
(higher harmonic contributions are neglected in this ap- 
proach). This result is model independent as long as 
the dominating interaction is J2 and the Fermi surfaces 
are located close to the F and M points in the folded 
BZ. The cosfcajCosfcj, form factor changes sign between 
the electron and hole pockets in the BZ. It resembles 
the order parameters of s* proposed from weak-coupling 
argumenta^^. 

The J2 coupling will be of particular importance in 
the following. We point out two key points on J2-related 
physics as it has appeared in the literature up to now. 
First, the effect of J2 is underestimated in most an- 
alytic models constructed based on the pure iron lat- 
tice with only onsite interactions since the J2 exchange 
coupling originates mostly from superexchange processes 
through As (P) or Sc (Te). Second, in the effective 
{t] — Ji — J2 model studied beforoii. the supcrconduct- 
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ing state is obtained only when the magnetic exchange 
couphng strength is of the same order as the hopping 
parameters (or the bandwidth) of the model. Therefore, 
as t > J, it requires the effective bands given by {i} be 
renormalized. However, the absence of double occupan- 
cies in the standard t — J model is not strictly imposed 
in such an intermediately coupled effective model where 
the bandwidth is assumed to be of similar order as the 
interaction scale. 

Comparing the predictions from weak coupling and 
strong coupling, the 122 iron chalcogenides provide an 
interesting opportunity to address the difference between 
the two perspectives. In this paper, we predict that 
the s-wave pairing symmetry is robust even in extremely 
electron-overdoped iron chalcogenides because the AFM 
J2 is the main factor for pairing and the Ji is ferromag- 
netic (FM), a conclusion drawn from both neutron scat- 
tering experiments^^— and the magnetic structure asso- 
ciated with 245 vacancy ordering^Si^. As we will show, 
the FM Ji significantly reduces the competitiveness of 
d-wavc pairing symmetry. We substantiate this claim by 
two different methods. First, wc solve the three orbital 
{t} — Ji — J2 model on the mean field level to show that 
the s-wave pairing is the leading instability regardless of 
the change of doping given that J2 is large. We calculate 
a full phase diagram as Ji varies from FM to AFM. If 
Ji is AFM, we obtain a SC state with a mixed s-wave 
and d-wave pairing. Second, we use the functional renor- 
malization group (FRG) to analyze this trend obtained 
by mean field analysis for a 5-band model of the chalco- 
genides. We confirm that a dominant AFM J2 generally 
leads to robust s-wave pairing while an AFM Ji tends 
to favor d-wave pairing in the electron overdoped region. 
The competition between s-wave and d-wave weakens the 
superconducting instability scale. In addition, it drives 
the anisotropy feature of the superconducting form fac- 
tor as consistently obtained for various weak coupling 
approaches. Together, our analysis provides an expla- 
nation for the different behavior of superconductivity in 
the iron pnictides and iron chalcogenides in the electron 
overdoped region since Ji has opposite signs for these 
two classes of materials, i. e. Ji is AFM in the iron 
pnictides^i^ and FM in the iron chalcogenides. Our 
study suggests that determining the pairing symmetry of 
the 122 iron chalcogenides can provide important insight 
regarding whether the local AFM exchange couplings are 
responsible for the high superconducting transition tem- 
peratures. 

The paper is organized as follows. In Section |lTl we 
present the mean field analysis of the t — Ji — J2 model 
to show the differences between the iron pnictide setup 
Ji > and the chalcogenide setup Ji < in the electron- 
overdoped regime. This is followed by FRG studies in 
Section IIIII where we mainly investigate the competition 
between s-wave and d-wave in the effective model, and 
also analyze the possible effect of an additional pocket 
at the r point of the unfolded BZ which wc find to fur- 
ther increase the robustness of the s-wave pairing. The 



qualitative trends confirm the results obtained in Sec. |TI1 
In Section IIVI we provide a combined view on the chalco- 
genides and point out that the ferromagnetic sign of Ji 
is important to explain the robustness of s-wave pair- 
ing symmetry in these compounds. Furthermore, we set 
our work into context of other approaches to the prob- 
lem. We conclude in Section |V] that electron-overdoped 
chalcogenides exhibit a robust s-wave pairing phase when 
the NNN interactions are correctly taken into considera- 
tion. 



II. MEAN FIELD ANALYSIS 

We calculate the mean-field diagram of an effective 
model for the AFe2Se2 compounds. As the main relevant 
orbital weight is given by the dxz, dyz, and dxy orbital, we 
employ a three-orbital kinetic model with Ji and J2 in- 
teractions. For the case of strong electron doping we are 
interested in, we do not find qualitative differences when 
four or five orbital models are used. For a more thorough 
discussion of these aspects, refer to Section HITl The spe- 
cific kinetic theory we use for the mean-field analysis is 
a modified three-band model^, given by 

(Tii{k)-^i Ti2(fc) T^m \ 
r(fc)= T2i(fc) T22{k)--fi , (1) 

V T3i(fc) T32(fc) n3{k)-flj 

where 

T'ii(fc) = 2t2 cos{kj;) + 2ti cos(fey) + 4^3 cos{kj;) cos{ky), 

T22{k) = 2ti cos{kj;) + 2t2 cos(fey) + At^ cos{kj;) cos{ky), 

Tssik) = 2tr,{cos{kx) + cos{ky)) + At^ cos[kx) cos{ky) + 5, 

Ti2(fc) = 4t4sin(/c,)sin(A:j,),T2i(fc) = T*^{k), 

Tiz[k) = 2it-! sin[kx) + Aits sm{kx) cos{ky), 

T23{k) ~ 2it-i sin(fcy) 4- 4*^8 sin(fcy) cos(A:a;). (2) 

The other matrix elements are given by hermiticity. 
The parameters in the model are taken to be t = 
(0.02,0.06,0.03,-0.01,0.35,0.3,-0.2,0.1), 6 = 0.4, and 
/i — 0.412. (Throughout the article, energies are given 
in units of eV unless stated otherwise). The parameter 
set chosen gives the Fermi surface shown in Fig. [T] with 
a filling factor of 4.41 electrons per site. Aside from a 
negligibly small electron pocket at the M point in the 
unfolded Brillouin zone, the main features are the large 
electron pockets at X which dictate the physics of the 
mean-field phase diagram at this electron doping regime 
(sec FigH]). In the three band model, the small electron 
pocket appears around the M-point in the unfolded Bril- 
louin zone which may be related to the resonance feature 
experimentally discussed for the F point in the folded 
zone. In contrast, the 5-band fit to the chalcogenides we 
employ in Section IIIII suggests small electron pocket fea- 
tures around the F-point of the unfolded Brillouin zone. 
Despite this discrepancy, later we will see that both ap- 
pearances have a similar effect and can hence be discussed 
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FIG. 1. (color online) The Fermi surface used to represent the 
chalcogenides. Colors indicate the orbital components: Red 
dxz, Green dyz, and B\ue=dxy. The A and B are auxiliary 
labels used in Fig. |6l 



tice, we can classify the pairing symmetries according 
to the one-dimensional irreducible representations of the 
C4 symmetry. Since the pairing is a spin singlet, we can 
classify them as follows: an order parameter is of A- type 
(B-type) if it is even (odd) under a 90-degree rotation. 
This classification leads to six candidate pairings with A- 
symmetry and another six candidates with B-symmctry 
as the SC pairings include NN from Ji and NNN from 
J2 bonds, which manifests as the A-typc symmetry 
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and the B-typc symmetry 



(4) 



on the same footing. The interaction part in our mean 
field analysis is the pairing energy obtained by decou- 
pling the magnetic exchange couplings^i, which can be 
written as 

~^J2bli^r,r+x+y^a,r,r+x+y + J2bl^,r^r+x-y^a,r,r+x-y) 



where 



\iCa,r'.'\- represent sin- 



Ay — ^1 < ba^r,r+y >, 
and A"_j^ J2 < ba,r,r 
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'-^x+y 



J2 < baj-^ 



r+x+y 



> 



+x+y >, where a denotes the 
orbital index and a;, y are the two unit lattice vectors. 
We only consider intra-orbital pairing and ignore inter- 
orbital pairing which is very small as shown in previous 
calculationaii. Considering the C4 symmetry of the lat- 
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(5) 



glct pairing operators between the r, r' sites. 

Before we perform the self-consistent mean field calcu- 
lation, we define the pairing order parameters as follows: 
in real space, the pairings on two NN bonds and two 
NNN bonds are represented by A^;! = Ji < ba,r,r+x >, where 



In reciprocal lattice space, the mean-field Hamiltonian is 
given by 



fr^^l f{k) A{k) 

^ \ A^k) ~f*{-k) J ' 



Aii(A:) 
A(fc) = I A22(fc) 

A33(A:) 



(6) 



(7) 



and 



Aii(fc) = (A^^_, + A^^_J(cos(/c,) + cos{ky)) + (A^^jv,d + A^jv,d)(cos(/c,) - cos(fcj,)) 

+ 2(A^^^ J. -f Afj^j^ g) cos(fca; ) cos(fcy) + 2{A^jy^^^ + Af ,^) sm{k.j;) sm{ky), 

A22{k) = {A^N.s - ANN.s)icos{k^) + cos{ky)) + (A^ ^ - A^^ ,i)(cos(A::,) - cos(A:y)) 

+2(A^^jv s — Af ,,) cos(fc2,) cos(fcy) + 2(Af ~ ^NNN.d) sin{kx) sm(ky), 

Assik) = A^''^_^(cos(fc^) + cos(fcy)) + A^^^__^(cos(A:a;) - cos(A:y)) + 2A'^j^^^^cos{k^} cos{ky} 

+2A^Vjv,<iSin(A::r)sin(fcy). (8) 

I 



We emphasize that in above definitions the symbols s and space, and do not correspond to whether pairing is even 
d merely represent the geometric factor of pairing in k- or odd under a 90-degree rotation in a multi-orbital sys- 
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FIG. 2. (color online) The mean field phase diagram along 
Ji = in the parameter space (up); the quasiparticle spec- 
trum at J2 = 1 (middle) and J2 = 0.75 (bottom). At J2 = 1, 
one can see that the quasiparticle spectrum explicitly breaks 
C4 symmetry because of the mixing of A- and B-type pairing 
symmetries. 



tern. In general, there are more than one self-consistent 
set of {A}'s as self-consistent meanfield solutions. The 
free energies in each solution hence have to be compared 
to find the solution with the lowest free energy. 

First consider pure NNN-pairings stemming from J2 
(this is a reasonable limit to start with since Ji in 
FeTe(Se) has been shown to be ferromagnetic, thus not 
contributing to pairing in the singlet pairing channel). 
J2 is increased from zero to J2 = 1.5 while the band 
width is ~ 4. The robust superconductivity solution 
with purely A-type s-wave pairing is obtained when J2 
is larger than 0.4. This is to say the pairing remains 
the same as in iron-pnictides with the geometric factor 
cos{kx) cos{ky)^. The Bogoliubov particle spectrum is 
completely gapped in this state. When J2 becomes larger 





FIG. 3. (color online) The mean field phase diagram along 
J2 = in the parameter space (up); the quasiparticle spec- 
trum at Ji = 0.75 (bottom). At Ji = 0.75 one can see that 
the quasiparticle spectrum explicitly breaks the C4 symmetry 
because of the mixing of A- and B-type pairing symmetries. 



than 1, the ground state is a mixture of A- and B-type 
pairings. The nonzero B-type pairings all have the geo- 
metric factor sin(fca;) sin(A:j,) (see the phase diagram show 
in Fig. [2]). In the coexistence phase, the quasiparticle 
spectrum shows nearly gapless features at several points, 
and moreover, the dispersion explicitly breaks C4 rota- 
tion symmetry (see Fig. [2] displaying the quasiparticle 
spectrum of the lowest branch) . 

Second, we study the phase diagram when only (anti- 
ferromagnetic) Ji is present. In this case, only NN pair- 
ings are nonzero and there arc six SC gaps. Wc increase 
Ji from Ji = to Ji = 1.5 where the band width is 
W ^ 4. The SC order becomes non-zero from Ji = 0.4 
on. However, in this case, the B-type SC order arises 
slightly earlier than A-type SC order. The ground state 
is always a mixture of A- and B-type pairings. The two 
leading orders are A-type s-wave and B-type d-wavc in 
xz, yz orbitals while the sub- leading ones are s- and d- 
waves in the xy orbitals (Fig. Due to strong mixing 
of A- and B-type pairings, the quasiparticle spectrum is 
very anisotropic. It is, however, still nodeless, in contrast 
to a pure s-wave pairing cos{kx) + cos{ky) where there 
are nodes^ on the electron pockets (see Fig. [3]). 

Finally, for Ji and J2 antiferromagnetic, we fix Ji -|- 
J2 = 1 and change Ji — J2 as a parameter. We observe 
that NNN pairings dominate for Ji — J2 < —0.1 and 
NN pairings dominate for Ji — J2 > 0.2 (Fig. 0]). In the 
intermediate range, there is only weak B-type pairing. A 



5 




FIG. 4. (color online) The mean field phase diagram with 
Ji + J2 = 1 in the parameter space. 




B N N s-wave + 
B NNN s-wave 



ANN s-wave J 



FIG. 5. (color online) A schematic phase diagram for the 
model ^ within < Ji, J2 < 1. 




FIG. 6. (color online) The orbital resolved pairing amplitude 
on the FS for a typical s-wave (d-wave) pairing state in the 
upper (lower) panel, calculated within meanfield approxima- 
tion. The interaction parameters are Ji — 0, J2 ~ 0.8 for 
the upper panel and Ji = 0.5, J2 = for the lower panel. 
In the left half of these figures, the k-point traces the elec- 
tron pocket around X-point counterclockwise from point A 
in Fig. [1] and in the right half, it traces the electron pocket 
around the Y-point counterclockwise from point B in Fig. [T] 



schematic phase diagram within the range < Ji , J2 < 1 
is shown in Fig. [S] 

In the whole parameter region of (Ji, J2), the SC or- 
der parameters always have the same sign for all three 
orbitals. This can be seen in Fig. [6] where the orbital 
resolved pairing amplitude is shown along electron pock- 
ets around X. This result is essentially consistent with 
the FRG result shown in Sec. Hn] (FiglH]). It is, how- 
ever, different from what one would expect from the very 
strong coupling limit: There, the strong inter-orbital re- 
pulsion favors different signs of pairing for the dxy orbital 
and the d^z/yz orbitals^. Some quantitative differences 
between Fig. [5] and Fig. [T^] may be explained as the in- 
completeness of a three-orbital model and the fact that 
the meanfield pairing is not constrained to the FS. In 
Fig. [5] we also see that the orbital resolved pairing am- 
plitude is highly anisotropic: This is a natural reflection 
of different orbital composition on different parts of the 



Fermi surface. 

Following the Fermi surface topology in Fig. [U these 
mean-field results have been obtained in the case where 
a small electron pocket was still present at the M-point. 
For completeness of the analysis, we also adapted the 
parameters such that the electron pocket around the M 
point vanishes, leaving two pockets around X. Without 
the M pocket, we see that s-wave pairings are less favored 
than before, as its geometric factor is cos(A:2;) cos(A:j,) 
or cosikx) + cos{ky), both being maximized around M. 
With the two-pocket FS, taking Ji = and increasing 
J2 > 0, B-type pairings blend in at smaller J2 than shown 
in Figl5] taking J2 = and increasing Ji > 0, A-type 
pairings appear at slightly larger Ji than shown in Fig. [S] 
The main features still remain unchanged. These trends 
are in accordance with the FRG studies in the following 
section. 
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III. FRG ANALYSIS 

To substantiate the mean field results above, we em- 
ploy functional renormalization group (FRG)^"— to fur- 
ther investigate the pairing symmetry of the t-Ji-J2 
model. As an unbiased resummation scheme of all chan- 
nels, the FRG has been extended and amply employed to 
the multi-band case of ir on p nictides. More details can 
be found in Refs. [l^ [13, l39l and [s^ The conventional 
starting point for the FRG are bare Hubbard-type inter- 
actions which develop different Fermi surface instabilities 
as higher momenta are integrated out when the cutoff of 
the theory flows to the Fermi surface. To address the spe- 
cial situation found in the chalcogenides where the Fe-Se 
coupling is strong, not only local, but also further neigh- 
bor interaction terms would have to be taken into ac- 
count: in our FRG setup, the onsite Hubbard-type inter- 
actions of the same type as used in the study of pnictides 
triggers no instability at reasonable critical scales. This 
suggests already at this stage that the chalcogenides may 
necessitate a perspective beyond pure weak coupling. In 
addition, the total parameter space of bare interactions is 
large and constrained RPA parameters are not yet avail- 
able for this class of materials. For the purpose of our 
study, we hence constrain ourselves to the i- Ji- J2 model 
from the outset. This implies that the pairing interaction 
is already attractive on the bare level, and a development 
of an SC instability is expected as the physics is domi- 
nated by the pairing channel. Still, we can employ FRG 
to investigate the properties and competition of different 
SC pairing symmetries for the chalcogenides for different 
(Ji, J2) regimes. 

Within FRG, we consider general J1-J2 interactions 
which are not limited to the spins in the same orbital: 

H = Ji"^ "^{Sia ■ Sjfc - -riianjb) 

{%,]) a,b 



> 
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FIG. 7. (color online) The spectrum and the Fermi surfaces of 
the band structure proposed by Maier et al.— , colored accord- 
ing to the dominant orbital content. The color code is red d^z, 
green dyz, blue d^cy, orange d^2_y2 and magenta d^^2_^2. The 
numbered crosses show the center of Fermi surface patches 
used in the FRG calculations. 



decompose the pairing channel into eigenmodes. 



(9) 



and obtain the band-resolved form factors of the leading 
and subleading SC eigenmode (i.e. largest two negative 
eigenvalues). This way we are able to discuss the inter- 
play of c?-wave and s-wave as well as the degree of form 
factor anisotropy for a given setting of (Ji, J2). Com- 
paring divergence scales Ac gives us the possibility to 
investigate the relative change of Tc as a function of 
{JiT^)- Furthermore, we also investigate the orbital- 
resolved pairing modes^ by decomposing the orbital four 
point vertex 
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ni ,. . . ,714 — 1 



|VA(fcl,ni;fc2,"^2;fc3,"'3;fc4,'^4) 
X (fclXn2 (fc2)Uc«3 (fc3)w<in4 (^4) | , (10) 

where the m's denote the different orbital components of 
the band vectors. By investigating the intraorbital SC 
pairing modes in (jlOp . we make contact to the findings 
from the previous mean field analysis. 



The kinetic theory will differ in the various cases studied 
below. For all cases, we will study the full 5-band model 
incorporating all Fe d orbitals. Concerning the discretiza- 
tion of the BZ, the RG calculations were performed with 
8 patches per pocket, and a lOradius x Sangic mesh on 
each patch. (This moderate resolution is convenient to 
scan wide ranges of the interaction parameter space; 
we checked that increasing the BZ resolution did not 
qualitatively change our findings.) The output of the RG 
calculation is the four-point vertex on the Fermi surfaces: 
VA(fci, ni; ^2,^2; fca, ng; ^4, "4)4^„,,43„35Cfc^„,,Cfc^„^5, 
where the fiow parameter is the IR cutoff A approaching 
the Fermi surface, and with ki to fc4 the incoming 
and outgoing momenta. We only find singlet pair- 
ing to be relevant for the scenarios studied by us: 
Efc,p^A(fc:P)[OlOp], where Of = c^.^C-fc^- We 



A. Two-pocket scenario 

We start by studying the 5-band model suggested be- 
fore by Maier et ali^i. There are only two electron pock- 
ets at the X point of the unfolded Brillouin zone closely 
resembling the Fermi surface topology and orbital con- 
tent employed for our mean-field analysis (Fig. [7|). 

The RG fiow and the form factors of the leading di- 
verging channels are shown in Fig. [5] for dominant J2 
and in Fig.[S]for dominant Ji. As stated before, the pair- 
ing interaction is already present at the bare level in the 
model so that we achieve comparably fast instabilities as 
the cutoff is flowing towards the Fermi surface. As found 
in Rcf . [Til , the dominant J2 scenario exhibits a leading s- 
wave cos cos ky form factor which causes the same sign 
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FIG. 8. (color online) Typical RG flows and the supercon- 
ducting gaps associated with the Fermi surfaces for the two- 
pocket scenario with (Ji, J2) — (0.1,0.5) eV. Leading form 
factor is denoted in blue, sub-leading form factor in green. 



FIG. 9. (color online) Typical RG flows and the supercon- 
ducting gaps associated with the Fermi surfaces for the two- 
pocket scenario with (Ji, J2) = (0.5,0.1) eV. Leading form 
factor is denoted in blue, sub-leading form factor in green. 



on both electron pockets (blue dots in Fig. |8]). The sub- 
leading form factor is found to be of d-wavc cos — cos ky 
type, changing sign from one electron pocket to the other. 
The inverse situation is found for dominant Ji . As shown 
in Fig.ini the d-wavc cos kx~cos ky form factor establishes 
the leading instability. As before, the form factor does 
not cross zero related to nodeless SC for this parameter 
setting. 

With only pairing information available on the limited 
number of sampling points along FS, it is impossible to 
obtain as in the mean-field analysis the superconducting 
gap in the whole BZ. For illustration, a mixture of a 
small A-type NNN d-wave pairing and large A- type NNN 
s-wave pairing is indistinguishable from a pure A-type 
NNN s-wave pairing; a mixed state of a small B-type NN 
s-wave pairing plus a large B-type NN d-wave pairing, 
and a state with pure B-type NN d-wave pairing show 
little difference if one compares the gap on a few points 
along the Fermi surfaces. For this reason, the symbol 

y2 used in this section refers to a pairing consisting of 
a large A-type NNN s-wave pairing and possible small 
components of A-type NNN d-wave pairing or A-type 
NN s/d- wave pairing. In turn, the symbol d^2_y2 refers 
to a pairing made up with a large B-type NN d-wave 
pairing and possible small components of B-type NN s- 
wave pairing or B-type NNN s/d- wave pairing. 

We have scanned a large range of (Ji, J2). For each 
setup we have obtained Ac as well as the ratio of the in- 
stability eigenvalues between s-wave and d-wave in the 
pairing channel (encoded by the two-color circles shown 
in Fig. [T0|) . The FRG result is qualitatively consistent 
with the mean field analysis. In the antiferromagnetic 
sector, the s-wave wins for dominant J2 while the d-wave 
wins wins for dominant Ji. For ferromagnetic Ji cor- 
responding to the situation in chalcogenides, we find a 
robust preference of s-wave pairing. The anisotropy of 
the s-wave gap around the pockets in the FRG calcu- 
lation is also qualitatively consistent with the meanfield 



result. The gap on the Fermi surfaces with dxy orbital 
character is smaller than the gap on those with dxz,yz 
orbital character. 

The predictions from the mean field analysis are fur- 
ther confirmed for the mixed phase regime where s-wave 
and d-wave coexist in the mean field solution. In FRG, 
one of these instabilities will always be slightly preferred; 
still, when both instabilities diverge in very close prox- 
imity to each other, this regime behaves similarly to the 
coexistence phase. For illustration, in Fig. [11] we have 
plotted the dependence of Ac on Ji — J2, with Ji + J2 
fixed to 0.7 eV; there is a clear reduction of the critical 
scale (and thus the transition temperature) when there 
is a strong competition between s- and d-wave channels. 

Following ((TU)) we also analyze the orbital decomposi- 
tion of the SC pairing from FRG (Fig. [T^. We constrain 
ourselves to the most relevant three orbitals d^y, dxz, and 
dyz- In particular, we observe that the SC orbital pair- 
ing induces the same sign for all three dominant orbital 
modes, in correspondence with the mean field analysis 
presented before. 



B. Three-pocket scenario 

Recent ARPES data^ on the chalcogenides may sug- 
gest the existence of a shallow flat pocket around the F 
point (the location, and especially the kz position of such 
a pocket are still under debate). By tuning parameters, 
we have obtained in the previous section a three-orbital 
model that has an additional electron-pocket around M- 
point in the unfolded BZ. 

In our FRG approach we can take a more profound mi- 
croscopic perspective on this issue. From the true band 
structure calculations at hand for the chalcogenides, we 
consider it unlikely that it will be a hole band regular- 
ized up towards the Fermi surface. Instead, we investi- 
gate the effect of a possible electron band at the F point 
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FIG. 10. (color online) The phase diagram of the two- 
pocket model. Each pie shows the relative strengths of the 
two leading pairing channels, with the radius proportional to 
[8 + logio(Ac/eV)]'^. The color code for pairing symmetries is 
green 8^,2^2 and red d^2_y2. 




FIG. 11. (color online) The variation of the critical scale 
Ac along a line through parameter space which interpolates 
between s and d wave. A minimum is visible for comparable 
ordering tendency in s-wave and d-wave. See the caption of 
Fig. [10] for more details on the pie charts. 
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FIG. 12. (color online) The orbital-resolved pairing form fac- 
tors of two typical RG flows. The upper resides in the dom- 
inant s-wave and the lower in the d-wave regime. The color 
code is the same as in Fig. [71 i. e. red d^z, green dyz, blue 
dxy, orange d^2_y2 and magenta d322_^2. 

was achieved by 

H^H+ C(k)cL«a(k)<(k)ckfc„ 

k,a.b,s 



in the unfolded Brillouin zone. This is suggested from 
the folded 10-band calculations, where one electron-type 
band closely approaches the Fermi level around the F 
poinli^. This band should be very flat and shallow. From 
the weak coupling perspective of particle-hole pairs cre- 
ated around the Fermi surfaces, this will probably have a 
small effect: particle- hole pairs will only be created up to 
energy scales of the depth of the electron band at the X 
point below the Fermi surface, providing some hole-type 
phase space for the electron band at F. In a ( Ji, J2) pic- 
ture, however, this may still significantly promote scat- 
tering along T -ir^ X , which may further stabilize the s- 
wave phase regime. We have hence developed a modified 
band structure designed for this scenario. There, we have 
bent down the band dominated by d^y in the two-pocket 
model band structure^ without changing its band vector 
and created an electron pocket around F, accordingly of 
mainly dxy orbital content (Fig. llSp . The band bending 



where u(k) is the eigenvector of the band dominated by 
dxy The shift of energy ^(k) was intentionally chosen 
such that the F pocket exhibits some nesting with the X 
electron pockets. 

The phase diagram is shown in Fig. [121 FRG results 
for typical scenario for the s-wave and d-wave regime are 
shown in Fig. [T3l and Fig. [Ml respectively. As suspected, 
the additional pocket strengthens the tendency to form 
an s-wave in the system, aside from exhibiting an addi- 
tional constant s-wave instability in a small regime for 
dominant Ji. 



IV. DISCUSSION 

The above calculations demonstrate that the s-wave 
pairing symmetry is always robust if the AFM NNN J2 is 
strong while a d-wave pairing can be strong if Ji is AFM 
for the electron overdoped region. Moreover, if both of 
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FIG. 15. (color online) The band structure and the Fermi 
surfaces of the modified band structure, colored according 
to the dominant orbital content. The color code is red d^z, 
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numbered crosses show the center of Fermi surface patches 
used in the FRG calculations. 



FIG. 13. (color online) Typical RG flows and the super- 
conducting gaps associated with Fermi surface for the three- 
pocket scenario with (Ji,J2) = (0.2,0.8) eV. Leading form 
factor is denoted in blue, sub-leading form factor in green. 
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FIG. 14. (color online) Typical RG flows and the super- 
conducting gaps associated with Fermi surface for the three- 
pocket scenario with (Ji, J2) = (0.9,0.3) eV. Leading form 
factor is denoted in blue, sub-leading form factor in green. 



them arc AFM, there is strong competition between the 
s-wave and d-wave pairing. When there are hole pock- 
ets, as shown beforeii, even in a range of Ji ^ J2, the 
contribution to pairing from Ji is much weaker than the 
one from J2. In that case, an AFM Ji will not gener- 
ate strong d-wave pairing so that the s-wave wins easily. 
From neutron scattering experiments, it has been shown 
that a major difference between iron-pnictidcs and iron- 
chalcogenides is that the NN coupling Ji change from 
AFM in the formeri^Si^ to FM in the latterii. In fact, 
Ji is rather strongly FM in the latter, which explains 
the high magnetic transition temperature (500 K) in the 
245 vacancy ordering state as shown in Ref . |5ll Combin- 
ing these results, we can partially answer the question 
regarding the different behaviors between iron-pnictidcs 
and iron-chacogenides in the electron-overdopcd region: 
why can the high SC transition temperature be achieved 
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FIG. 16. (color online) The phase diagram of the three-pocket 
model. Here we are able to resolve the s-wave channel into 
constant s-wave (grey), extended s^2y2-wave (green), and the 
nodal s^2^y2-wave (purple). Parameter sets with Ji ~ — 1 
and J2 ~ have highly oscillating form factors which are due 
to artifacts in the calculation; the triplet channel would have 
to be considered in these cases. 



in the latter, but not in the former? Since Ji in iron- 
pnictidcs is AFM while it is FM in iron-chacogenides, Ji 
will weaken the SC pairing in the former but not in the 
latter. 

A few remarks regarding this work follow: (i) Our 
mean-field result is qualitatively consistent with the re- 
sults from a similar model with five orbitals^"— . The 
critical difference is regarding Ji being FM, and has not 
been addressed previously; (ii) s-wave pairing symmetry 
has also been obtained in Refs. IH, Q, and [g^- How- 
ever, the s-wave pairing only shows up either in a narrow 
region or with drastically different parameter settings. 
Therefore, the s-wave is not robust from a microscopic 
point of view. Instead, the c?-wavc is a robust result 
in these studies. Still, even the c?-wave pairing strength 
based on the scattering between two electron pockets is 
generally weak, as specifically discussed in^, which is 
another difficulty for this type of mechanism, (iii) Our 
results suggest that there is no difference between iron- 
pnictides and iron-chalcogenides in terms of pairing sym- 
metry. Both of them are dominated by s-wave pairing. 
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If both hole and electron pockets are present, the signs 
of the SC order in hole and electron pockets are oppo- 
site, namely s*. However, the mechanism causing 
is different in the weak and strong coupling approach. 
In the weak coupling approach, the sign change is due 
to the scattering between the hole and electron pockets 
while in the strong coupling approach, the sign change is 
due to the form factor of the SC order parameters which 
is specified to be cos kx cos ky since the pairing mainly 
originates from the AFM J2. Therefore, to obtain 
pairing symmetry, the existence of both hole and elec- 
tron pockets is necessary in the weak coupling approach, 
but not in the strong coupling one. (iv) The reason that 
the superconductivity vanishes in the iron pnictides in 
the electron overdoped region is not solely due to the 
competition between s-wavc and d-wave pairing symme- 
try. It is also due to the weakening of local magnetic 
exchange coupling themselves and the reduction of band 
width renormalization. (v) The prospective experimental 
confirmation of s-wave pairing symmetry in KFe2Se2 will 
support that superconductivity in iron-based supercon- 
ductors might be explained by local AFM exchange cou- 
plings, (vi) Neutron scattering also suggests that there 
is significant AFM exchange coupling between two third 
nearest neighbor sites, i.e. J^^^. The existence of J3 
will further enhance the s-wave pairing since it gener- 
ates pairing form factors as cos 2kx -\- cos 2ky in recipro- 
cal space which in turn can enhance the pairing at the 



electron pockets. 

V. CONCLUSION 

In summary, we have shown that the pairing symmetry 
in electron-overdoped iron-chalcogenides is a robust s- 
wave. The fact that the NN magnetic exchange coupling 
is FM, which diminishes the possibility of d-wave pairing 
symmetry in these materials. From a unified perspective 
of high- Tc cuprates and high-Tc chalcogcnides, the NN 
AFM exchange coupling gives rise to the robust d-wave 
pairing in the cuprates while the NNN AFM exchange 
coupling gives rise to the robust s-wave pairing in iron- 
based chalcogenide superconductors. 
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